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ABSTRACT 


Molecular dynamics simulation was performed to investigate the effect of thickness on the properties and struc- 
ture of ultrathin tetrahedral amorphous carbon films (ta-C). The present simulation showed that both the density 
and residual compressive stress decreased with increasing the thickness of ultrathin ta-C films, which was in 
agreement with the experimental results. The gradient of properties with the thickness was dependent on the 
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incident kinetic energy of deposited atoms; when the incident kinetic energy was 10 eV/atom, the slower 
gradient was observed. Further structural analysis indicated that the critical relaxation of highly distorted bond 
angles was responsible for the reduction of residual compressive stress in the films deposited at 10 eV/atom, 
while the joint relaxation of both the distorted bond angles and C- C bond length led to the fast release of residual 
compressive stress in the films deposited at 70 eV/atom. 


1. Introduction 


Owning to the excellent properties such as high hardness, low 
coefficient of friction, wear-resistance, extremely smoothness and 
chemical inertness, diamond-like carbon (DLC) films have been used 
as protective coatings on the disks and read-write heads of magnetic 
disk storage devices to minimize the mechanical wear and corrosion 
[1-3]. Recently, the storage density of hard disk drives is increasing at 
a very rapid rate. In order to meet the requirement for the high storage 
density, the magnetic spacing which is the vertical distance between the 
read head and magnetic storage layer, must be reduced. One effective 
way to reduce the magnetic spacing is to decrease the thickness of 
protective films [4,5]. So a protective layer with several nanometers is 
required to satisfy the demand of storage density increasing. Generally, 
ultrathin tetrahedral amorphous carbon (ta-C) films with unique com- 
bination of desirable mechanical properties and tribological lubricity 
are studied [5-7]. For example, Casiraghi et al. prepared ta-C films by fil- 
tered cathodic vacuum arc technology and showed that the film with 
the thickness of 2 nm was pin-hole free and had a Young's modulus of 
~100 GPa, an sp? content of ~50% and a roughness of ~0.12 [5]. Beghi 
et al. also found that ultrathin 2-nm-thick carbon-based films retained 
a Young's modulus of about 100 GPa [6]. 
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However, the narrow range of applications of ultrathin ta-C films, 
mainly in the fields of magnetic recording media, restricts the wide, in- 
depth studies on ultrathin ta-C films. In addition, it is known that the 
structure and properties of DLC films strongly depend on the thickness 
of films [8-10]. Especially for ultrathin ta-C films, the decrease of thickness 
leads to the complexity in experimental characterization and the strong 
relationship between the structure and properties [11]. Thus, in order to 
understand the evolution of structural properties with the thickness and 
provide guidelines for the preparation of the ultrathin ta-C films used in 
magnetic storage devices, a further insight into the ultrathin ta-C films 
from the viewpoint of atomic scale is required. Molecular dynamics simu- 
lation (MD) technique provides a robust method for the fundamental un- 
derstanding of the atomic bond structure and properties of DLC films in 
atomic scale. By MD simulation, previous studies have investigated the 
variations of structure and properties with the kinetic energy of deposited 
carbon atoms [12], the energy of the carrier gas ions such as Ar” [13] or 
growth species [14]. But the thickness dependence of the properties and 
structure of the ultrathin ta-C films is still not fully understood yet. 

In the present work, we deposited the ultrathin ta-C films with the 
thickness of 2-6 nm by a classical MD to investigate the effect of thick- 
ness on properties and structure. The dependence of residual stress and 
density on thickness of films was studied, and the structural evolutions 
including bond angles and bond lengths were mainly analyzed. Struc- 
tural analysis revealed that with the increase of thickness of ultrathin 
ta-C films, the distorted atomic bond structure was relaxed, which 
accounted for the reduction of residual stress. However, the incident 
energy of deposited carbon atoms also played a key role on the change 
of structural properties with the thickness of ultrathin ta-C films. 
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2. Computational methods 


Classical MD simulation was used to study the deposition process of 
ultrathin ta-C films with different thickness using the three-body em- 
pirical potential Tersoff [15]. Previous results have revealed the 
limitation of Tersoff potential where the m bonding is not adequately 
considered [16-18], and also developed Brenner [19,20], the 2nd gener- 
ation reactive empirical bond order (REBO) potential [21] and the adap- 
tive intermolecular reactive empirical bond order (AIREBO) potential 
[22] to simulate amorphous carbon growth, but Tersoff is still proved 
to be an effective and accurate potential for carbon-based systems. 

A diamond (001) single crystal was served as substrate with the 
dimensions of 25.2210 x 25.2210 x 24.0758 A? in the x, y and z direc- 
tions, which consisted of 28 atomic layers with 100 carbon atoms per 
layer and was equilibrated at 300 K for 100 ps before deposition. The 
incident carbon atoms were introduced at the position of 10 nm above 
the substrate surface at a random {x, y} position. The positions of atoms 
in the bottom two monolayers were frozen to mimic the bulk substrate, 
while all the other atoms were unconstrained. The deposition was sim- 
ulated using NVE ensemble implemented in the large-scale atomic/ 
molecular massively parallel simulator (LAMMPS) code [23]. Periodic 
boundary conditions were applied along the lateral xy dimensions. 

The kinetic energy of incident carbon atoms was fixed at 70 eV/atom, 
because it was the optimum energy for a highly stressed and dense ta-C 
film deposition [24]. For comparison, the kinetic energy of 10 eV/atom 
was also considered. The thickness variations of ultrathin ta-C films 
were achieved by changing the deposition time. The time step of 1 fs 
was used. The time interval between two sequential deposited carbon 
atoms was 10 ps. The previous report has indicated that the time inter- 
val of 10 ps was enough for relaxing the atomic structure and 
diminishing the unrealistic effect of high carbon flux on the deposition 
process [24]. However, a number of processes such as diffusion or rear- 
rangement processes that occurred on much longer time scales were 
neglected in the simulation due to the required computation time. The 
substrate temperature was rescaled to 300 K by the Berendsen method 
[25] after the atomic rearrangement caused by the bombardment of in- 
cident atoms was finished, and the heat bath coupling constant was 
10 fs in the simulations. 


3. Results and discussion 


The final configurations of the deposited films as a function of the 
thickness are summarized in Fig. 1, in which the red numbers corre- 
spond to the thickness of films calculated by subtracting the original 
thickness of substrate from the total thickness of systems after 
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Fig. 2. RDF of the deposited films with the different thickness when the incident kinetic 
energy is (a) 10 eV/atom or (b) 70 eV/atom, respectively. 


deposition. It shows that a typical amorphous film is generated for 
each case, which will be described later by the radial distribution func- 
tions (RDF). As the deposition time increases from 20 to 50 ns, the thick- 
ness of films deposited at 10 eV/atom varies from 23.19 to 57.97 A 
(Fig. 1a), while the case at 70 eV/atom increases from 24.08 to 57.97 A 
(Fig. 1b). It is deduced that changing the incident kinetic energy has lit- 
tle effect on the deposition rate of ultrathin ta-C films; Fig. 1 also shows 
that the surface roughness is independent on the thickness. In addition, 
it should be noted that when the incident kinetic energy is 70 eV/atom, 
the incident atoms can easily implant into the diamond lattice, resulting 
in an obvious intermixing layer between the film and substrate (Fig. 1b). 

The RDF, g(r), is proportional to the density of atoms at a distance r 
from another atom, which is a very important parameter for structural 
characterization of amorphous materials, and it is computed by the 
formula 


dN 


s(n) = p- 4nr?dr 


(1) 


where p is the average density of system, dN is the number of atoms 
from r to r + dr. Fig. 2 shows the RDF for all ultrathin ta-C films, in 
which the red dotted lines represent the positions of the 1st nearest 
neighbor (1.54 A) and 2nd nearest neighbor (2.52 A) of crystalline 
diamond, respectively. For all films with the different thickness, the 


57.97A 


Fig. 1. Cross-sectional snapshots of the deposited films with the different thickness when the incident kinetic energy is (a) 10 eV/atom or (b) 70 eV/atom, respectively. The red numbers on 
top of each film correspond to the thickness of films which was deposited by changing the time from 20 to 30, 40, and 50 ns, the colors from blue to red represent the height of the atoms, 


and the red dotted line means the position of substrate surface. 
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Fig. 3. Variations of density and hybridization ratio (sp*%, sp?% or sp%) along the growth direction at the thickness of (a) 46.37 A in the film deposited at 10 eV/atom and (b) 45.48 A in the 
film deposited at 70 eV/atom. The positive direction in horizontal axis corresponds to the films, while the negative direction represents the original substrate. 


typical amorphous character is observed, that is long-range disorder 

and short-range order. Increasing the thickness of ultrathin films 

makes no obvious effect on the positions of the 1st and 2nd peaks 
which are located at approximately 1.48 A and 2.53 A, respectively. 

Fig. 3 shows the variations of density and hybridization ratio (sp*%, 

> sp?% or sp%) along the growth direction of films deposited at 10 eV/atom 

and 70 eV/atom. Similar to previous study [24], the whole system is also 

© divided into four regions including substrate, transition region, stable 

region and surface region. It can be noted that the sp? fraction in all 

LO films appears low values, which might attribute to the limit of Tersoff 

potential where the n bonding was not adequately considered. Because 

© the density and hybridization ratio exhibit relative constant values in 

the stable region, the representative properties (density and compres- 

sive stress) as a function of the thickness of ultrathin ta-C films are 

quantified using this region, as illustrated in Fig. 4. The residual stress, 

o, is computed by the formula 


0 


(2) 


s= where d is the dimensionality of the system (2 or 3 for 2d/3d), V is the 
system volume (or area in 2d); and Pxx, Pyy and P,, are the diagonal com- 

C ponents of the stress tensor. It shows that both the density and residual 
compressive stress decrease with increasing the film thickness whatev- 

s= er the kinetic energy of incident carbon atoms. Notably, for the ultrathin 
= ta-C films deposited at 70 eV/atom, the density and compressive stress 
©) with the thickness decrease linearly; as the thickness increases from 
24.08 to 57.97 A, the compressive stress drops about 63.2%, while the 
density is only reduced by 8.5%. However, at the incident kinetic energy 

of 10 eV/atom, the density is linearly reduced from 2.83 to 2.63 g/cm? as 

the thickness varies from 23.19 to 57.97 A, but the gradient is slower 
than that of the films deposited at 70 eV/atom; the compressive stress 
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decreases from 5.47 to 1.79 GPa and the gradient with the thickness in- 
creases gradually. 

For comparison and validation with the simulation, the experiments 
were also carried out by preparing the ultrathin ta-C films with different 
thickness using a home-made 45° double-bent filtered cathodic vacuum 
arc (FCVA) technique. P-type silicon (100) wafers were used as sub- 
strates. Graphite with purity of 99.999% was used as the cathodic target. 
During deposition, the films were deposited at a bias voltage of — 80 V 
and the arc current was fixed at 60 A. The deposition time was changed 
from 4 to 20 min. The experimental results of ultrathin ta-C films with 
the thickness ranging from 7.6 to 24 nm are shown in Fig. 5. It reveals 
that following the thickness of ultrathin ta-C films the compressive 
stress also decreases significantly with a little deterioration of density. 
This is consistent with the simulation results. 

In order to elucidate the evolution of properties with the thickness of 
ultrathin ta-C films, the atomic bond structure needs to be further ana- 
lyzed. The distributions of both the bond angles and bond lengths for all 
ultrathin ta-C films are plotted to gain the fractions of distorted bond 
angles and bond lengths, as shown in Fig. 6. It is known that the stable 
bond angle and bond length of graphite are 120° and 1.42 A, respective- 
ly (red dotted line in Fig. 6); for crystalline diamond, the stable bond 
angle and bond length are 109.471° and 1.54 A (black dotted line in 
Fig. 6). Previous study has revealed that the high residual compressive 
stress is mainly originated from the distortion of both the bond angles 
and bond lengths of carbon network, which are less than 109.471° 
and 1.42 A, respectively [24]. By integrating the bond angle and bond 
length distributions (Fig. 6), the fractions of distorted bond angles 
(<109.471°) and bond lengths (<1.42 A) in each film are thus deduced 
separately, as shown in Fig. 7. 

First, it reveals that the films deposited at 70 eV/atom have higher 
fractions of both the distorted bond angles and bond lengths than that 
at 10 eV/atom, which results in the higher residual stress shown in 
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Fig. 4. Evolutions of (a) Density and (b) Compressive stress as a function of thickness of the films deposited at the incident kinetic energies of 10 eV/atom and 70 eV/atom, respectively. 
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Fig. 5. Density and compressive stress as a function of thickness of the films deposited by 
FCVA method. The residual stress was calculated from the curvature of the film/substrate 
composite using Stoney's equation and the XRR was conducted to characterize the density 
of films. 


Fig. 4b. For the films deposited at 10 eV/atom (Fig. 7), increasing the 
thickness from 23.19 to 34.78 A leads to the opposite changes in the 
fractions of distorted bond angles and bond lengths, and as a conse- 
quence which leads to the mediate change in the residual stress 
(Fig. 4b). Nevertheless, as the thickness increases from 34.78 to 
57.97 A, the fractions of distorted bond angles and bond lengths 
decrease simultaneously, and this combined contribution induces the 
fast reduction of residual stress, but the relaxation of distorted bond an- 
gles is dominated in this case. However, when the incident kinetic 
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energy is 70 eV/atom, Fig. 7 shows that the fractions of both the 
distorted bond angles and bond lengths as a function of the thickness 
decrease gradually, which accounts for the significant release of the 
residual stress (Fig. 4b). 

For the films deposited at 10 or 70 eV/atom, the different behaviors 
of atomic bond structure with the thickness can be understood by the 
impact of carbon atoms with high incident energy. When the kinetic en- 
ergy is 70 eV/atom, more energy with increasing the thickness of films 
is released to the atoms around the deposition area, which not only re- 
laxes the bond angles, but also effectively relaxes the C- C bond length 
with high strength and directionality (Fig. 7). Therefore, in the ultrathin 
ta-C films deposited at 70 eV/atom, the reduction of compressive stress 
is derived from the relaxation of both the highly distorted bond angles 
and bond lengths, which gives explanation for the significant gradient 
of residual compressive stress with the thickness in the simulation 
and experiments. 


4. Conclusion 


The ultrathin ta-C films with the thickness ranging from 2 to 6 nm 
were deposited by MD simulation using Tersoff potential. The effect of 
thickness on properties and structure was evaluated. For the films de- 
posited at 10 eV/atom or 70 eV/atom, the similar deposition rate was 
observed. The properties of ultrathin ta-C films seriously depended on 
the thickness. As the thickness of films deposited at 70 eV/atom varied 
from 24.08 to 57.97 A, the compressive stress decreased linearly from 
12.1 to 4.45 GPa with only 8.5% discount of the density. With the inci- 
dent kinetic energy decreased to 10 eV/atom, the minimal compressive 
stress of 1.79 GPa was obtained with the density of 2.63 g/cm? when the 
thickness was 57.97 A. This behavior agreed well with the experimental 
results of ta-C films deposited by FCVA technique. By analyzing both the 
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Fig. 6. Distributions of bond angle and bond length as a function of the thickness of films deposited at the incident kinetic energies of (a) 10 eV/atom and (b) 70 eV/atom, respectively. 
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Fig. 7. Fractions of (a) distorted bond angles (<109.471°) and (b) bond lengths (<1.42 A) as a function of the thickness of films deposited at the incident kinetic energies of 10 eV/atom and 


70 eV/atom, respectively. 


bond angle and bond length distributions, it revealed that when the in- 

cident kinetic energy was 10 eV/atom, the reduction of residual com- 

pressive stress with the thickness mainly attributed to the relaxation 

of the distorted bond angles, while in the films deposited at 70 eV/atom, 
Y™ it resulted from the relaxation of both the distorted bond angles and 

bond lengths. The study for the thickness dependence of properties of 
Y ultrathin ta-C films provided a theoretical guide for designing and 
©) preparing ultrathin ta-C films with high performance. 
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